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具有 学 习 机 制 的 正弦 余弦 算法 


方 坦 阳 ， 武 相 军 '， 游 大 涛 
(河南 大 学 软件 学 院 , 河南 开封 475004) 


摘 要 : 针对 标准 正弦 余弦 算法 在 求解 函数 优化 问题 时 易 陷入 局 部 最 优 、 收 敛 精度 较 差 等 问题 ， ose 
习 机 制 的 正弦 余弦 算法 。 该 算法 引入 精英 反 向 学 习 策略 构造 精英 及 反 向 群体 ， 1 ee 择优 保留 ， 从 而 优 
化 了 种 群 中 的 个 体位 置 、 提 高 了 算法 的 寻 优 精度 ; 同时 ， 利 用 个 体 的 反思 学 习 能 力 来 防止 个 体育 目地 向 当前 最 优 解 
学 习 ， 使 算法 停滞 在 局 部 最 优 ， 从 而 有 效 地 避免 了 算法 的 未 成 熟 收 化 。 通 过 13 个 标准 测试 菠 数 进行 仿真 实验 ， 实 
验 结果 证 明 ， 该 算法 相 比 于 对 比 算 法 具有 较 强 的 彰 棒 性 和 函数 优化 能 
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Sine cosine algorithm with learning mechanism 


Fang Xuyang, Wu Xiangjun', You Datao 
(College of software, Henan University, Kaifeng Henan 475004, China) 


Abstract: In order to solve function optimization problems that standard sine cosine algorithm is easy to fall into local 
optimum and poor convergence accuracy, the paper proposed a sine cosine algorithm with learning mechanism. The 
algorithm introduced the elite opposite learning strategy to construct the elite and the opposite population, and preserve the 
optimal individual in the mixed population, thus optimize the individual position in the population and improve the 
optimization accuracy of the algorithm. At the same time, the algorithm used individual reflective learning ability to avoid 
individuals blindly learning from the current optimal solution, and made the algorithm stagnating in the local optimal, thus 
effectively avoid the immature convergence of the algorithm. The algorithm performed simulation experiments with 13 
standard test functions. The experimental results showed that the algorithm was more robust and better function 
optimization ability than comparison algorithms. 
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0 ”外 j 言 子 计 算 与 正弦 余弦 算法 相 结 合 的 方法 ， 利 用 量子 比特 对 位 置 

进行 编码 、 量 子 旋转 门 对 个 体 状 态 进行 更 新 、 量 子 非 门 实现 

在 科学 和 工程 领域 ， 优 化 问题 一 直 备 受 关注 。 由 于 问题 变异 操作 等 ， 实 验 验证 了 算法 的 有 效 性 ， 达 到 较 好 的 效果 。 

的 复杂 性 ， 使 用 传统 方法 解决 有 时 是 不 可 能 前 ， 许 多 文献 [8] 引 入 反 向 学 习 的 方法 ， 对 当前 个 体 产 生 反 向 解 ， 扩 大 

研究 人 员 一 直 致 力 于 开发 一 种 新 的 优化 算法 ， 来 解决 现实 中 了 对 解 空 间 的 探索 。 文献 [9] 提 出 了 混合 灰 狼 优 化 的 正弦 余弦 

复杂 的 优化 问题 。 其 中 群体 智能 算法 应 用 较为 广泛 ， 如 粒子 算法 ， 利 用 正弦 余弦 更 新 公式 改进 头 狼 的 移动 方向 和 速度 ， 

群 优化 (particle swarm optimization，PSO) LI、 人 工蜂 群 算 使 算法 得 到 较 好 的 提升 。 文 献 [10] 提 出 了 结合 差分 进化 的 正 

法 (artificial bee colony algorithm，ABC) 外、 引力 搜索 算法 ” 弦 余 弦 算 法 ， 使 用 差分 进化 算 子 帮助 算法 跳出 局 部 最 优 解 区 
(gravitational search algorithm，GSA) BI 和 鲸鱼 优化 算法 域 。 

(whale optimization algorithm，WOA) 中 等 。 鉴于 此 ， 为 了 提高 SCA 算法 的 收敛 速度 和 避免 早熟 现 


正弦 余弦 算法 (sine cosine algorithm, SCA ) 1 是 Mirjalili 
于 2016 年 提出 的 一 种 新 型 的 群体 智能 优化 算法 ,该 算法 结构 


象 , 本 文 在 基本 的 SCA 算法 上 提出 了 精英 反 向 学 习 策略 与 反 


思 学 习 策略 相 结合 的 改进 方法 。 通 过 对 标准 测试 函数 的 仿 


简单 、 参 数 较 少 且 易 于 实现 ， 它 的 搜索 过 程 主要 受 正弦 和 余 ”实验 ， 并 与 基本 的 SCA、PSO、WOA 及 其 几 种 改进 的 SCA 
弦 函 数 的 影响 。 然 而 ， 像 其 他 群体 智能 优化 算法 一 样 ， 存 在 。 ”算法 相 比 较 ， 实 验 表 明 本 文 所 改进 的 算法 显著 提高 了 算法 的 
收敛 精度 低 、 收 敛 速度 慢 、 易 陷入 局 部 最 优等 问题 。 近 两 年 。 全 局 优化 能 力 。 

和 许多 研究 人 员 从 不 同 的 角度 提出 了 不 同 的 改进 方法 ， 以 PN 

提高 算法 的 优化 能 力 。 文 献 [6] 提 出 了 一 种 转换 参数 非 线 性 递 1 正弦 余弦 算法 

减 的 正弦 余弦 算法 ， 分 别 以 抛物 线 函 数 和 指数 函数 控制 参数 SCA 算法 是 最 近 提 出 的 一 种 新 颖 的 全 局 优化 算法 ， 有 别 
i 的 变化 ， 实 验 结果 表明 指数 函数 对 参数 i 的 调整 ， 能 够 较 ”于 根据 不 同 生物 启发 机 理 的 群体 智能 优化 算法 。 在 SCA 算法 
好 的 平衡 算法 的 全 局 勘探 和 局 部 开发 能 力 。 文献 [7] 提 出 了 量 。 中 , 个 体 的 更 新 状态 主要 依赖 于 两 个 数学 函数 Sine 和 Cosine 
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的 变化 来 实现 优化 搜索 。 假设 第 t 代 种 群 中 个 体 X! 由 D 个 分 
量 组 成 ， 即 Xf =(xhy x Kh)，ie{l,2..., N}, de{l,2,..., D}, 
N 为 种 群 规模 ，D 为 个 体 维度 。 种 群 的 初始 化 使 用 式 (1) 随 
机 产生 个 体 所 在 位 置 。 
Xr! =X4,, +rand(0,1)( Xu — Xu,,) (1) 

其 中 X% 和 Xs 表示 个 体 在 第 d 维度 的 上 下 界 。 
在 每 次 迭代 中 ， 第 i 个 个 体 的 位 置 ， 将 按 以 下 更 新 方程 
移动 所 在 位 置 ， 即 

Xi+nxsin(n )xInBs —X!| nn<0.5 
> 下 + oo n>0.5 


(2) 
Bw 表示 种 群 的 


其 中 :表示 在 第 t 次 欠 代 中 个 体 ; 的 位 置 ; 


当前 最 优 位 置 ; 是 服从 均匀 分 布 的 随机 数 ， 有 e[0,27]， 
ss[-223]，2s[04，5 是 控制 参数 。 


如 式 2) 所 示 ， 式 中 有 四 个 主要 的 调节 参数 分 别 是 、 
请、 和 nn。 参 数 决定 下 一 次 迭代 时 第 i 个 个 体 的 位 置 区 


域 , 定义 了 下 次 过 代 的 最 大 可 能 移动 的 区 域 范围 ，» 定义 了 
当前 解 应 该 朝 目 标 解 前 进 还 是 远离 ， 参 数 是 目标 解 的 一 个 
随机 权 值 ， 当 |sl>1 增 强 目 标 解 对 当前 解 的 影响 ， 当 |s|<1 目 


标 解 对 当前 解 的 影响 减弱 ; ” 是 正弦 和 余弦 机 制 切换 的 随机 
概率 。 Msin(5) 或 icos(5) 共同 领导 着 算法 进行 局 部 搜索 和 全 
局 搜索 ， 当 Msin(5) 或 5cos(5) 的 值 大 于 1 或 小 于 -1 时 ， 进 行 
全 局 勘探 ， 在 -1 与 1 之 间 时 ， 进 行 局 部 开发 。 为 了 使 算法 的 
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个 体 在 搜索 空间 中 存在 一 个 反 向 个 体 ， 反 向 个 体 可 能 更 接近 
于 最 优 位 置 。 这 样 就 能 在 每 次 搜索 过 程 中 ， 同 时 搜索 当前 个 
体 和 当前 个 体 的 反 向 解 ， 对 其 进行 择优 保留 。 保 留 下 来 的 优 


秀 个 体 加 速 了 种 群 中 个 体 向 最 优 位 置 靠 近 ， 提 高 了 算法 的 求 
解 效 率 。 
定义 1 反 向 解 (opposite solution)H2。 设 在 区 间 [lb,ub] 


上 存在 一 个 实数 x， 则 x 的 反 向 数 定义 为 x=Ib+ub-x。 鉴于 
此 ， 假 设 在 R 域 上 存在 一 个 点 X 在 D 为 的 空间 上 ， 则 
且 ss ， 即 它 的 反 向 解 
和 =( 交 ap) ,由 式 必 =k*(15+ub) 玉 生成 ,为 [0,1] 之 间 均 勾 
分 布 随 机 数 ， 称 为 一 般 化 系数 。 

定义 2 基于 反问 解 的 优化 (opposite-based 
optimization)03。 设 存在 种 群 P， 其 反 向 种 群 为 P' ,基于 优 者 
保留 的 方式 从 PuP' 的 种 群 中 选择 最 优秀 的 n 个 个 体 体 组 成 
新 种 群 ， 称 为 基于 反 向 解 的 优化 。 

定义 3 动态 一 般 反 向 学 习 (dynamic generalized 
opposition-based learning)021。 设 存在 种 群 X， 数 量 为 N， 维 
度 为 D， 则 第 t 次 迭代 过 程 中 ， 种 群 的 反 向 解 为 
Xi =k*(1b; tuby)—Xi， 其 中 ib;=min(X;) , ub;=max(X;) 。 Xi 为 


种 群 第 i 个 体 在 第 j 维 上 的 分 量 。 
定义 4 ”精英 反 向 解 (elite opposition solution)02。 设 在 D 


X=(N,%,..., Xp) ， 


勘探 和 开发 能 力 得 到 平衡 ， 寻 找到 解 空间 中 可 能 的 区 域 ， 
最 终 收 敛 到 全 局 最 优 ，' 通过 式 (3) 进行 自 适应 调整 ， 即 


a 
1 (3) 


其 中 : + 是 当前 迭代 的 次 数 ，7 是 最 大 迭代 次 数 ;，a 是 一 个 
常数 。 在 算法 中 ，4a 一般 设 置 为 2。 

SCA 算法 通过 数学 函数 的 变化 来 改变 初始 状态 的 位 置 ， 
种 群 中 的 个 体 更 新 依靠 函数 值 的 增加 或 者 减少 来 随机 地 更 新 
每 次 迭代 中 每 个 个 体 的 状态 ， 使 种 群 在 前 期 保持 多 样 性 ， 后 
期 随 着 ;的 减少 ， 个 体 趋 于 局 部 开发 ， 最 终 收敛 于 全 局 最 优 
状态 。 算 法 的 基本 原理 如 图 1 所 示 。 


© 


{ 


图 1 正弦 余弦 算法 的 基本 原理 


Fig.1 Basic principle of sine-cosine algorithm 


2 ”改进 的 正弦 余弦 算法 


2.1 精英 反 向 学 习 

SCA 算法 在 迭代 后 期 将 对 当前 全 局 最 优 位 置 附 近 进 行 
很 小 的 邻 域 内 搜索 ， 不 断 尝 试 更 新 最 优 解 。 如 果 搜 索 过 程 远 
离 了 理论 最 优 解 ， 算 法 就 很 难 在 短 时 间 内 收敛 到 全 局 最 优 。 
对 此 本 文 引入 精英 反 向 学 习 策 略 ,扩大 精英 种 群 的 搜索 范围 ， 
从 而 改善 了 SCA 算法 收敛 速度 慢 、 精 度 不 足 等 问题 。 
2.1.1 基本 概念 

反 向 学 习 (opposition-based learning) 0 22 策略 是 近年 
来 智能 优化 邻 域 出 现 的 一 种 新 技术 ， 目 前 已 在 人 工蜂 群 算法 
(ABC)、 差 分 进化 算法 (DE)、 粒 子 群 算法 PSO) 等 群体 


维 空 间 上 Xu=(02 sz) 为 当前 种 群 的 精英 个 体 
Xi =(00pop) 的 反 向 解 ， 该 反 向 解 定 义 为 
有 ;=k*(1b+ub)-x ,ke[0,1] 为 服从 均匀 分 布 的 随机 数 ， 利 用 该 
系数 可 以 生成 精英 个 体 的 多 个 反问 解 。 
通过 生成 的 精英 反 向 解 ， 增 加 了 种 群 向 全 局 最 优 收敛 的 
有 益 信 息 ， 可 以 加 强 对 最 优 个 体 周围 邻 域 的 探索 ， 提 高 算法 
的 局 部 开发 能 
2.1.2 精英 反 向 策略 
为 了 减少 算法 偏离 全 局 最 优 位 置 的 可 能 性 ， 加 强 对 优秀 
个 体 周 围 空间 的 搜索 是 非常 有 必要 的 ， 这 种 改进 可 以 提高 算 
法 勘探 新 解 的 能 力 。 本 文 将 精英 反 向 学 习 的 策略 融入 SCA 算 
法 中 ， 充 分 利用 精英 种 群 的 信息 ， 搜 索 精英 个 体 及 其 反 向 解 
所 在 的 空间 。 具 体操 作 如 下 : 
种 群 中 的 个 体 由 式 (2) 执行 后 ， 对 种 群 进行 排序 ， 取 其 
中 10% 的 优秀 个 体 组 成 精英 种 群 B 。 
计算 个 体 XePis 的 边界 [oab ,并 求解 其 动态 边界 
[min(b;), max (ub;)] 。 
根据 定义 3、4 生成 个 体 Xs 的 动态 精英 反 向 种 群 已， 。 
车 反 向 种 群 Bs 超出 了 Xs 的 动态 边界 [min(1b;),max(ub;)] 
的 限制 ， 则 由 式 (1) 在 边界 内 随机 生成 新 个 体 进行 替换 。 
将 混合 种 群 BUR 根据 适应 度 大 小 从 高 到 低 排 序 ， 选 
择 前 1/2 的 个 体 进 入 下 一 代 。 
循环 执行 bp~e, 直 到 达到 停止 条 件 ， 算 法 结束 。 
2.2 反思 学 习 策 略 
在 教 与 学 优化 算法 [3 与 中 , 所 有 的 学 生 通 过 教师 的 课堂 
教学 学 习 科 目 知识 ， 同 时 还 向 身边 同学 请 教 来 不 断 提 高 自己 
的 学 习 能 力 。 文 献 [16] 提 出 基于 自主 学 习 行为 的 教 与 学 优化 
算法 ， 把 学 习 分 类 为 自主 学 习 阶段 和 反思 阶段 。 由 于 学 生 个 
体 无 法 充分 发 挥 自身 的 学 习 能 力 ,最 终 影响 整 体 的 学 习 效 果 ， 
所 以 在 求解 复杂 函数 优化 问题 时 ， 容 易 过 早 陷 入 局 部 最 优 区 
域 。 受 此 启发 ， 本 文 对 SCA 算法 进行 以 下 改进 ， 提升 
学 习 能 力 。 


在 SCA 算法 中 , 种 群 的 个 体 只 依赖 于 当前 最 优 解 来 进行 


智能 优化 算法 中 得 到 了 成 功 的 应 用 。 它 的 主要 思想 是 : 每 个 


自身 状态 的 更 新 。 因 此 ， 算 法 有 较 大 的 可 能 陷入 局 部 最 优 状 
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态 ， 导 致 了 算法 无 法 寻找 到 满意 的 解 。 此 时 ， 需 要 对 个 体 进 
行 局 部 变异 操作 ， 个 体 以 当前 的 学 习 结果 向 前 一 次 的 学 习 情 
况 进 行 反 思 学 习 ， 增 加 逃逸 局 部 区 域 的 概率 。 反 思 学 习 的 公 
式 如 (4) 所 示 ， 即 


Xi = Xi+OO( Xi —X!) (4) 
其 中 : Xi 表示 在 第 t 次 迭代 中 个 体 i 的 位 置 ，Xi 表示 执行 式 
(2) 后 的 位 六 表示 经 过 反思 过 程 产生 的 新 位 置 ，@% 表 
示 学 习 因 子 ，2s[-11 ;， @ 表示 点 乘 。 
每 一 代 个 体 在 执行 正弦 或 余弦 操作 后 ， 开 始 进入 自我 反 
思 阶 段 。 根 据 自身 前 一 次 的 学 习 状 态 反思 当前 的 学 习 状 态 ， 
对 当前 学 习 状 态 进行 自我 修正 。 通 过 自身 学 习 的 方式 能 够 更 
好 的 防止 个 体 易 受 当前 全 局 最 优 解 影响 ， 而 陷入 局 部 最 优 。 
这 里 为 了 防止 反思 过 程 中 的 随机 性 太 大 , 将 用 式 (5) 对 学 习 
妹子 。 进行 约束 ， 同 时 为 了 避免 学 习 能 力 退 化 、 增 强 算 法 的 
改 敛 性 ， 使 用 贪 禁 学 习 的 方式 ， 个 体 在 反思 前 后 的 学 习 状 态 
根据 学 习 成 绩 《〈 适 应 度 ) 进行 择优 选择 。 
W= C7) xcos(n) (5) 
其 中 : 是 [0,*] 上 的 随机 数 ; C 是 常数 ， 经 实验 C=100 时 效果 


小 ” 
MY 
习 


子 随 迭代 次 数 变化 情况 如 图 2 所 示 。 


0 200 400 600 800 1000 
和 迭代 次 数 

图 2 ”学习 因子 随 迭 代 次 数 变化 情况 
Fig.2 Learning factor changes with number of iterations 
图 2 可 以 看 出 ，C0" 在 迭代 过 程 中 非 线 性 递减 ， 即 个 
体 在 前 期 具有 较 强 的 反思 活动 ， 实 现 自我 的 超越 ， 提 高 勘探 
能 力 ; 在 后 期 随 着 adenal 趋 于 全 局 最 优 状态 
反思 活动 不 断 减少 ， 学 习 因 子 。 的 强度 逐渐 变 弱 。 
在 式 (6) 中 ， 个 体 在 通 过 反思 后 ， 若 优 于 反思 前 的 个 体 
将 对 该 个 体 进行 蔡 换 。 使 用 贪 禁 选 择 的 方式 将 最 优 个 体 送 入 
下 一 代 ， 确 保 了 每 一 代 个 体 都 趋向 于 全 局 最 优 解 方向 移动 。 

el 而 


X; Otherwise 


其 中 : f() 表示 个 体 的 适应 度 。 

反思 学 习 策 略 能 够 较 大 地 提高 算法 的 全 局 优化 能 力 ， 个 
体 对 自己 当前 的 学 习 状 态 进行 深入 思考 ， 并 根据 上 次 学 习 状 
态 进 行 变 异性 调整 ， 有 效 地 防止 了 自己 停滞 在 局 部 状态 ， 提 
高 了 规避 局 部 最 优 的 能 
2.3 改进 的 算法 步骤 

为 了 便于 后 续 的 实验 对 比分 析 改 进 后 的 算法 记 为 
MSCA(modified sine cosine algorithm)。 改 进 后 的 算法 以 伪 代 
码 来 表示 ， 算 法 的 具体 流程 如 表 1 所 示 。 

MSCA 算法 : 

1: 初始 化 参数 ， 利 用 式 〈1 ) 随机 产生 种 群 数量 为 N 的 个 体 

和 一 ( XXXp); 

2: 计算 每 个 个 体 的 适应 度 ; 
3: 设置 最 优 解 的 位 置 及 适应 度 ; 
4: do 
0 


for i=ltoN 


更 新 让 万, 三 ,六 
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gE 证 <0.5 
8: Xi =X!+nxsin(n)x|nB, — XH 
9: else 

10: Xi =X!+nxcos(n)x|nDos —XH| 


11: 正弦 余弦 操作 后 产生 的 个 体 的 新 位 置 进 行 反 思 学 习 
X= Xr +OO(X} —X!) 
12: 选择 反思 前 后 适应 度 最 好 的 个 体 作为 下 一 代 : 
wil 
本 


Otherwise 


L3: end for 
14: 对 当前 种 群 的 前 10% 的 最 优 个 体 计算 反 向 位 置 : 


X; =rand*(lb; +ub; )—X: 


15: 选择 1/2 的 Bs Bs 混合 群体 中 的 最 优 个 体 ; 
16: 计算 种 群 适应 度 ,更 新 最 优 解 位 置 及 适应 度 ; 


17: while (current iteration < maximum iterations) 


18: return best Solution obtained as the global optimum 

MSCA 算法 主要 由 初始 化 、 正 弦 余 弦 位 置 更 新 、 反 思 学 
习 、 精 英 反 向 学 习 四 部 分 组 成 。 由 上 述 算法 的 流程 可 得 ， 假 
设 算法 的 迭代 次 数 为 m， 种 群 数 量 为 n, 维 度 为 4， 其 中 种 群 
初始 化 和 个 体 的 位 置 更 新 与 基本 的 SCA 算法 复杂 度 相同 ,为 
0O(mnd) 。 在 反思 学 习 过 程 中 ， 对 每 个 个 体 进行 位 置 扰动 ， 其 
时 间 复 杂 度 为 0(mna) ; 在 精英 反 向 学 习 的 过 程 中 ， 需 要 对 
10% 的 个 体 进 行 反 向 计算 且 构 造 混合 群体 进行 比较 ， 其 时 间 
复杂 度 为 0(mnd) 。 因 此 ， 整 个 算法 的 时 间 复 杂 度 为 0(mnd) 。 


3 ”实验 仿真 


3.1 实验 环境 及 测试 函数 

为 了 验证 MSCA 算法 对 函数 的 优化 能 力 , 本 文采 用 文献 
[9] 中 的 13 个 标准 测试 函数 进行 实验 对 比分 析 。 仿 真实 验 运 
行 环境 为 Intel Core i7 CPU， 主 频 3.60 GHz， 内 存 8 GB， 
Windows7 64 位 操作 系统 ， 实 验 仿真 软件 采用 MATLAB 
R2017b。 实 验 中 采用 的 算法 统一 设置 为 种 群 规 模 N=30， 最 
大 和 迭代 次 数 T=1000， 算 法 在 每 个 标准 测试 函数 上 独立 运行 
30 次 ,统计 其 平均 值 和 标准 差 。 对 比 的 算法 和 描述 如 表 1 所 
示 。 测 试 函数 如 表 2、3 所 示 。 

表 1 对 比 算法 及 描述 


Table 1 Contrastalgorithms and descriptions 
算法 描述 
SCA 标准 正弦 余弦 算法 
OBSCAB 基于 反 向 学 习 的 正弦 余弦 算法 
HGWOSCA®™! 混合 灰 狼 优 化 的 正弦 余弦 算法 
WOA 标准 鲸鱼 优化 算法 
PSO 标准 粒子 群 算法 
里 中 . PSO 的 参数 设 为 C=2, C=2, Wa =0.9, W,, =0.2 。 
表 2 单 峰 基准 测试 函数 


Table 2 Unimodal benchmark functions 


Fun 函数 名 搜索 区 间 论 最 优 值 
Fl Sphere [-100,100] 0 
F2 Schwefel 2.22 [-10,10] 0 
F3 Schwefel 1.2 [-100,100] 0 
F4 Schwefel 2.21 [-100,100] 0 
F5 Rosenbrock [-30,30] 0 
F6 Step [-100,100] 0 
F7 Quartic Noise [-1.28,1.28] 0 
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表 3 多 峰 基准 测试 函数 (b) ~ (d) F10、F12 和 Fl13 函数 的 收敛 结果 可 以 看 出 ， 该 
Table 3 Multimodal benchmark functions 函数 易 陷 入 局 部 最 优 , 但 MSCA 最 终 的 寻 优 精 度 较 其 他 算法 
Fun 函数 名 搜索 区 间 理论 最 优 值 仍然 表现 较 好 。 
Fs Schwefel 2.26 [-500,500] -418.9829*n 表 4 函数 FL~EF7 的 寻 优 精度 对 比 
F9 Rastrigin [-5.12,5.12] 0 ‘Table 4 Variance and mean information of six algorithms( Fl~F7) 
请 0 Ackley [-32,32] 0 on 30 维 100 维 
Fll Griewank [-600,600] 0 平均 值 标准 差 ” 排序 ”平均 值 标准 差 ”排序 
F12 Penalized1l [-50,50] 0 WOA 1.3535E-150 4.7989E-150 2 7.9419E-147 4.3008E-146 2 
Fl3 Penalized2 [-50,50] 0 OBSCA 7.0799E-28 2.6100E-27 4 8.0438E-09 2.6023E-08 4 
3.2 实验 结果 分 析 HGWOSCA 1.3165E-50 1.8525E-50 3 1.6487E-25 1.4927BE-25 3 
为 了 体现 本 文 算 法 在 问题 维 数 和 复杂 度 上 的 优越 性 ， 实 下 SCA 6.2969E-03 1.0909E-02 6 4.7299E+03 3.5584E+03 6 
验 将 分 别 对 维 数 为 30、100 两 种 不 同 的 函数 维度 进行 测试 。 PSO 8.4698E-09 2.0700E-08 5 2.5551E+00 1.4886E+00 5 
其 中 FL~EF7 为 连续 性 单 峰 值 函 数 ，F8~FI3 为 连续 性 多 峰 函 MSCA ~ 1.2395E-233 0.0000E+00 1 2.5194E-220 0.0000E+00 1 
数 ， 均 是 具有 普遍 性 、 代 表 性 的 常用 于 检测 群体 智能 算法 全 WOA 2.7529E-103 1.0638E-102 2 7.0937E-103 3.1981E-102 2 
局 优化 能 力 的 函数 。 函 数 的 表达 式 可 以 查阅 文献 [9]。 OBSCA 5.1910E-26 1.6873E-25 4 2.4920E-13 6.7044E-13 4 
1) 寻 优 精度 分 析 HGWOSCA 6.9748E-30 1.0040E-29 3 1.2850E-15 5.0447E-16 3 
表 4 列 出 了 六 种 算法 在 单 峰 值 函数 求解 时 30 次 运行 结果 gcA 16508B.05 285436.05 5 27108E100 4.2881E100 5 
的 平均 值 和 标准 差 的 对 比 。 由 表 4 可 见 ，MSCA 算法 在 Fl 函 PSO 3.5512E-04 8.2176E-04 6 1.1212E+01 4.1496E+00 6 
数 上 的 求解 精度 高 出 WOA 算法 及 其 他 改进 的 正弦 余弦 算法 MSCA 9.7071E-108 5.3011E-107 1 2.4938E-106 1.3045E-105 1 
80 ~200 个 数量 级 ， 对 F3、F4 函数 的 求解 精度 较 对 比 算法 效 WOA 23101E+04 1.0649E+04 6 9.0586E+05 1.7474E+05 6 
果 十 分 显著 ; 在 F6 函数 的 求解 上 ，MSCA 算法 在 高 维 的 求 OBSCA 2.5799E-03 1.2901E-02 3 7.3762E+03 6.8831E+03 3 
解 能 力 优 于 较 好 的 PSO 算法 ， 在 低 维 的 求解 上 较 PSO 算法 HGWOSCA 1.1660E-12 3.7569E-12 2 1.9520E+01 4.2611E+01 2 
结果 差 ， 在 F5、F7 函数 上 ， 算 法 对 30、100 维 的 求解 精度 1 
较 低 ， 但 MSCA 算法 较 其 他 五 种 算法 仍 具 有 优势 ; 故 总 体 考 PSO 1.6050E+01 8.5760E+00 4 1.0830E+04 2.4641E+03 4 
虑 ，MSCA 在 单 峰 函数 问题 的 求解 精度 、 算 法 稳定 性 (标准 MSCA = 9.0933E-151 4.8622E-150 1 4.3847E-114 2.2422E-113 1 
差 ) 等 方面 比 对 比 算法 具有 更 好 的 效果 。 WOA 3.1947E+01 2.8662E+01 6 7.0368E+01 2.7823E+01 5 
表 5 列 出 了 六 种 算法 在 多 峰值 函数 求解 时 30 次 运行 结果 OBSCA 7.4303E-05 1.5470E-04 3 3.8343E+01 1.1390E+01 4 
的 平均 值 和 标准 差 的 对 比 。 由 表 5 可 以 看 出 , 对 于 F8~Fl3 多 HGWOSCA 1.0879E-11 2.3705E-11 2 1.7681E-03 2.2282E-03 2 
峰值 函数 的 寻 优 结果 , 在 问题 维度 为 30、100 维 时 MSCA 均 We gCA 20936E101 1.1576B101 5 85688B101 3.3002E100 6 
优 于 对 比 算法 ;在 F8 函数 的 求解 上 , 相 比 于 对 比 算法 , MSCA PSO 6.2625E-01 1.6545E-01 4 8.3142E+00 1.1723E+00 3 
能 够 找到 接近 理论 最 优 解 的 结果 ,SCA 在 该 函数 上 的 处 理 能 MSCA 1.7361E-113 9.0691E-113 1 9.8756E-112 5.3552E-111 1 
力 最 差 ， 在 F9、F11 函数 的 求解 上 ，WOA、MSCA 算法 能 WOA 2.7015B+01 4.0893E-01 2 9.7582E+01 4.3771E-01 2 
达到 理论 最 优 值 且 MSCA 在 30、100 维 求解 性 能 较为 稳定 ; OBSCA 2.8197E+01 2.5030E-01 4 9.9191E+01 7.8873E-01 4 
在 F10、F12、F13 函数 的 求解 上 ，MSCA 的 寻 优 结果 与 理论 HGWOSCA 2.7151E+01 9.1918E-01 3 9.7905E+01 6.2451E-01 3 
最 优 解 相 比 还 存在 精度 不 足 ， 但 优 于 对 比 算法 。 可 见 MSCA gcA 73237E102 29100B103 6 53408E107 23688E107 6 
在 多 峰值 函数 上 的 优化 能 力 比较 强 ， 逃 逸 局 部 最 小 值 的 能 PSO 4.9549E+01 3.1024E+01 5 2.4044E+03 1.3275E+03 5 
相 比 而 言 要 好 于 对 比 算法 。 故 总 体 考虑 ，MSCA 对 多 峰 函 数 MSCA 1.9731E-05 2.9621E-05 1 5.5978E-05 9.1330E-05 1 
问题 的 求解 表现 具有 很 大 的 优越 性 。 WOA 5.7319E-02 6.8223E-02 3 2.0318E+00 7.9950E-01 2 
2) 收敛 曲线 分 析 OBSCA 43876E+00 1.8999E-01 5 2.2373E+01 4.6430E-01 5 
为 了 进一步 考查 算法 的 平均 进化 趋势 , 实验 中 列 出 了 30 HGWOSCA 1.1532E+00 4.4087E-01 4 1.0407E+01 1.1248E+00 4 
维 时 算法 的 收敛 曲线 。 限 于 篇 幅 限 制 ， 在 文中 仅 列 出 了 F1、 gcA 471208100 4.41806-01 6 69219B103 39950B403 6 
F5、F6、F7 单 峰 测试 函数 ，F9、F10、F12、F13 多 峰 测 试 函 PSO 5.2557E-09 8.2397E-09 1 2.7054E+00 1.2978E+00 3 
数 的 平均 函数 值 的 迭代 收敛 曲线 。 部 分 单 峰值 函数 的 收敛 曲 MSCA 5.0650E-08 9.6047E-08 2 2.1809E-07 3.5146E-07 1 
线 如 图 3 所 示 。 由 图 3 (a) ~ (d) 可 以 看 出 ，MSCA 算法 的 WOA 2.1976E-03 2.5725E-03 3 1.5681E-03 1.5729E-03 2 
收敛 曲线 下 降 最 快 ， 并 能 够 随 着 迭代 次 数 的 增加 持续 收 钱 ， OBSCA 2.4692E-03 1.5650E-03 4 1.0853E-02 6.2737E-03 4 
出 现 停滞 现象 。 其 中 WOA 算法 是 较 新 的 优化 算法 ， 它 在 F7 HGWOSCA 1.0902E-03 5.9821E-04 2 3.6505E-03 1.3202E-03 3 
单 峰 值 函 数 的 求解 精度 相 比 其 他 群体 智能 算法 要 好 。 通 过 图 SCA = 2.8870E-02 2.3546E-02 5 6.6319E+01 4.0981E+01 5 
中 显示 , MSCA 算法 在 Fl1、F5、F7 函数 上 收敛 趋势 优 于 WOA PSO 6.7384E-02 2.7413E-02 6 1.4099E+03 3.6454E+02 6 
算法 ， 也 优 于 其 他 改进 的 SCA 算法 。 从 图 3 (b) 可 知 ， 对 MSCA 3.5682E-05 3.2475E-05 1 4.3021E-05 3.8397E-05 1 
比 算法 在 后 期 收敛 出 现 停 滞 ， 收 敛 精度 较 低 ,而 MSCA 算法 以 上 收敛 图 表明 了 MSCA 算法 可 以 以 更 快 的 收敛 速度 
能 规避 停滞 发 生 。 由 图 3 (c) 可 看 出 ，WSCA 与 PSO 相 比 得 到 精度 更 高 的 优化 结果 ,在 六 种 算法 中 全 局 搜索 能 力 最 佳 。 
在 后 期 迭代 效果 不 佳 ， 前 期 收敛 下 降 较 快 ， 但 优 于 图 中 的 对 3) 综合 实验 分 析 
比 算法 。 仿真 实验 的 算法 排序 结果 如 表 6 所 示 。 表 中 列 出 了 六 种 
图 4(a) ~(d) 是 部 分 多 峰值 函数 的 收敛 曲线 。 从 MSCA ”算法 在 13 个 标准 测试 函数 的 30、100 维度 的 平均 排序 结果 及 


算法 的 收敛 性 上 可 以 看 出 ， 在 F9 函数 上 收敛 速度 明显 ,在 。” 总 排序 。 在 同等 的 实验 前 提 下 ，MSCA 算法 相 比 于 对 比 算法 
F10、F12 和 F13 函数 上 ， 也 比 对 比 算法 收敛 程度 好 。 由 图 4 
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有 显著 优势 。 
表 5 函数 F8~Fl3 的 寻 优 精度 对 比 
Table 5 Variance and mean information of six algorithms( F8~ Fl3 ) 
30 维 100 维 

Mn 全 二 平均 值 “标准 差 ”得 平均 值 标准 差 站 
厅 订 

WOA -11549.179 1.2447E+03 2 -3.7602E+044.5384E+03 2 
OBSCA -4168.9017 3.1396E+02 5 -7.4608E+034.9664E+02 5 
HGWOSCA -5827.5 7.4939E+02 4 -1.6342E+041.7700E+03 4 
四 SCA -3998.299 3.3629E+02 6 -7.2570E+034.6787E+02 6 
PSO -6451.3732 1.0114E+03 3 -1.8731E+046.1457E+03 3 
MSCA -12569.483 5.7187E-03 1 -4.1898E+041.1390E-02 1 
WOA 0.0000E+000.0000E+00 1 0.0000E+00 0.0000E+00 1 
OBSCA 0.0000E+000.0000E+00 1 2.7440E-06 1.4944E-05 2 
HGWOSCAS.3355E+007.8261E+00 2 1.1795E+01 8.2861E+00 3 
a SCA 1.0882E+011.5187E+01 3 2.3894E+02 1.0306E+02 4 
PSO 4.9957E+011.0100E+01 4 4.5198E+02 5.4017E+01 5 
MSCA 0.0000E+000.0000E+00 1 0.0000E+000.0000E+00 1 
WOA 3.9672E-15 2.7572E-15 2 3.3751E-15 2.3137E-15 2 
OBSCA 2.1618E-02 9.0653E-02 5 4.7178E+00 3.7659E+00 5 
HGWOSCA 2.7060E-14 3.6712E-15 3 2.7930E-13 3.3242E-14 3 
四 SCA 1.5319E+018.0031E+00 6 1.9520E+01 3.3554E+00 6 
PSO 6.1330E-05 7.4154E-05 4 2.6409E+00 3.2792E-01 4 
MSCA 8.8818E-160.0000E+00 1 8.8818E-16 0.0000E+00 1 
WOA 2.3940E-03 1.3113E-02 3 0.0000E+00 0.0000E+00 1 
OBSCA 3.2839E-13 1.7985E-12 2 4.0223E-08 1.4392E-07 2 
HGWOSCA 2.7258E-03 5.3645E-03 4 5.4231E-03 1.1424E-02 3 
SCA 2.2328E-01 2.8470E-01 6 4.0373E+01 2.8646E+01 5 
PSO 8.6210E-03 1.0422E-02 5 4.8571E-02 3.0360E-02 4 
MSCA 0.0000E+000.0000E+00 1 0.0000E+000.0000E+00 1 
WOA 7.5227E-03 5.9677E-03 3 1.7488E-02 6.5294E-03 2 
OBSCA 4.9258E-01 6.6286E-02 5 1.1934E+00 1.8625E-01 4 
HGWOSCA7.6204E-02 4.2358E-02 4 3.2058E-01 7.7470E-02 3 
4 SCA 2.5894E+021.0303E+03 6 1.6412E+08 1.0791E+08 6 
PSO 6.9113E-03 2.6302E-02 2 2.2695E+00 1.1890E+00 5 
MSCA 4.7934E-097.7293E-09 1 1.9856E-09 3.8542E-09 1 
WOA 2.5202E-01 2.0052E-01 3 1.5783E+00 6.7054E-01 2 
OBSCA 2.3644E+001.6501E-01 5 1.1341E+01 8.3623E-01 4 
HGWOSCA 9.2799E-01 3.1218E-01 4 6.8008E+00 4.4165E-01 3 
3 SCA 1.1109E+012.4752E+01 6 3.3122E+08 1.8469E+08 6 
PSO 2.5637E-03 4.7266E-03 2 1.3012E+017.5262E+00 5 
MSCA 9.4655E-08 1.8124E-07 1 1.1467E-07 1.4137E-07 1 
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了 算法 寻找 最 优 解 的 机 会 实验 结果 表明 ，MSCA 算法 
的 改进 策略 是 有 效 的 且 具 有 较 强 的 鲁 棒 性 和 全 局 优化 能 力 ， 
有 望 应 用 于 复杂 问题 的 求解 。 同 时 ，SCA 作为 一 种 较 新 的 群 
体 智能 优化 算法 具有 很 大 的 工程 应 用 空间 ， 有 待 今后 进一步 
的 研究 。 
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表 6 算法 排序 结 
Table6 Algorithm sorting result 


算法 SCA OBSCA HGWOSCA WOA PSO MSCA 


平均 排序 5.50 3.85 3.04 2.65 4.19 1.04 
总 排序 ”6 4 3 2 5 1 
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